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■ On the basis of the percolation picture from the hadronic phase with hyperons to 



<N 



the quark phase with strangeness, we construct a new equation of state (EOS) with the 



O ■ pressure interpolated as a function of the baryon density. The maximum mass of neutron 



stars can exceed 2Mq if the following two conditions are satisfied; (i) the crossover from 
the hadronic matter to the quark matter takes place at around three times the normal 
nuclear matter density, and (ii) the quark matter is strongly interacting in the crossover 
CO ■ region. This is in contrast to the conventional approach assuming the first order phase 

transition in which the EOS becomes always soft due to the presence of the quark matter 
at high density. Although the choice of the hadronic EOS does not affect the above 
conclusion on the maximum mass, the three-body force among nucleons and hyperons 
^ plays an essential role for the onset of the hyperon mixing and the cooling of neutron 

^ . stars. 

Subject Index Neutron stars, Nuclear matter aspects in nuclear astrophysics, Hadrons and 
quarks in nuclear matter. Quark matter 
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OO . 1. Introduction 
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Neutron star (NS) is a cosmic laboratory which provides us with a testing ground for the 
rich phase structure of quantum chromodynamics (QCD) [1] through the observables such as 
^ I the mass (M), the radius (i?), the surface temperature (Tg), the surface magnetic field (Bg) 

and so on [2]. Among others, M and R are particularly important probes for constraining 
the equation of state (EOS) and the composition of high density matter. 

From the theoretical point of view, the onset of the strangeness degrees of freedom inside 
the NSs has attracted much attention in recent years: General consensus is that the hyperons 
(y) such as A and would participate in NS cores at densities of several times nuclear 
matter density (po = 0.17 fm-^)[3-9]. The precise value of the threshold density pth depends 
on the hyperon-nucleon interactions which have still uncertainties at the moment but will 
be improved by the future hypernuclear data [10-12] and by the lattice QCD simulations 
[13]. From the observational point of view, a massive NS, PSR J1614-2230, with Mobs — 
(1.97 ± O.O4)M0 was recently discovered [14]. Conflict between the 2M0-NS which requires 
stiff EOS and the y-mixing which gives soft EOS leads to a challenging problem whether 
massive neutron stars are in contradiction to the existence of the exotic components such as 
the hyperons and deconfined quarks[15-23]. 



'''These authors contributed equally to this work. 



© The Author(s) 2012. Published by Oxford University Press on behalf of the Physical Society of Japan. 
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 
(http://creativecommons.0rg/licenses/by-nc/3.O), which permits unrestricted use, 
distribution, and reproduction in any medium, provided the original work is properly cited. 



The purpose of the present paper is to investigate whether the "hybrid stars" which have 
quark matter in the core are compatible with 2M0-NS. Historicahy, the transition from the 
hadronic matter to the quark matter has been assumed to be the first-order phase transition 
and the Gibbs phase equihbrium conditions are imposed. However, treating the point-hke 
hadron as an independent degree of freedom is not fuhy justified in the transition region 
because ah hadrons are extended objects composed of quarks and gluons. Furthermore, the 
system must be strongly interacting in the transition region, so that it cannot be described 
neither by an extrapolation of the hadronic EOS from the low-density side nor by an extrap- 
olation of the quark EOS from the high-density side [24]. This is analogous to the BEC-BCS 
crossover realized in the many-body system of ultra-cold fermionic atoms [25]. 

Fig. 1 illustrates the above situation in terms of the pressure as a function of baryon density 
(p). One may expect a gradual onset of quark degrees of freedom in dense matter associated 
with the percolation of finite size hadrons, i.e., a smooth crossover from the hadronic matter 
to the quark matter. Such a percolation picture of hadrons has been discussed in seminal 
works such as Refs.[26, 27]. Also, hadron-quark continuity [28, 29] and hadron-quark crossover 
[30, 31] have been discussed in relation to the existence of color superconductivity at high 
density. As we will show in this paper, the crossover picture does not necessarily lead to the 
softening of EOS unlike the case of first-order transition, so that it provides us with a novel 
mechanism to support massive neutron stars with quark core. Preliminary account of our 
results has been reported in [32]. 

This paper is organized as follows. In §2, the characteristic features of the hadronic EOSs 
(H-EOSs) to be used at low densities are summarized. In §3, we treat the strongly interact- 
ing quark matter by using the Nambu-Jona-Lasinio (NJL) type model and derive the quark 
EOS (Q-EOS) to be used at high densities. In §4, we describe our interpolation procedure 




Fig. 1 Schematic picture of the QCD pressure (P) as a function of the baron density (p) under the assumption 
of the hadron-quark crossover. The crossover region where finite-size hadrons start to overlap and percolate is shown 
by the shaded area. The pressure calculated on the basis of the point-like hadrons (shown by the dashed line at low 
density) and that calculated on the basis of weakly interacting quarks (shown by the dashed line at high density) lose 
their validity in the crossover region, so that the naive use of the Gibbs conditions by extrapolating the dashed lines 
is not justified in general. 
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to obtain the EOS in the hadron-quark crossover region. In §5, numerical results and discus- 
sions are given for the bulk properties of hybrid stars, such as the M — R relationship, the 
maximum mass Mmax and the M — pc (central density) relationship. We discuss how these 
results depend on the different choice of H-EOS and Q-EOS. A comment on the cooling 
of NSs with respect to the hyperon mixture inside the core is also given. §6 is devoted to 
concluding remarks. 

2. Hadronic EOS (H-EOS) 

We consider several different EOSs with F-mixing: 

o TNI2, TNI3, TNI2u and TNISu [33, 34]: TNI2 and TNI3 are obtained by the G-matrix 
calculation with Reid soft-core potential for NN and Nijmegen type-D hard-core poten- 
tial for YN and YY. Also, a phenomenological three-body force [35] is introduced in a 
form of effective NN force to reproduce the saturation point of symmetric nuclear matter 
with the incompressibility A^=250MeV (TNI2) and /^=300MeV (TNI3). For TNI2u and 
TNI3u, the three-body interaction is introduced universally in a form of effective NN^ 
NY and YY forces. 

o AV18+TBF and Paris+TBF [36]: They are obtained by the G-matrix calculation but 
with different choice of potentials, AV18 and Paris potentials for NN and Nijmegen soft- 
core potential for YN and YY. Also, the three-body force of Urbana-type is introduced 
in a form of effective NN force to meet the saturation condition. 

o SCL3AI1 [39]: It is based on a relativistic mean field (RMF) model with chiral SU(3) sym- 
metry and logarithmic potential motivated by the strong coupling lattice QCD approach. 
Phenomenological parameters of the model are determined to reproduce the satura- 
tion condition, bulk properties of normal nuclei and separation energies of single- and 
double-A hypernuclei. 

In Fig. 2, we plot the pressure P for the F-mixed neutron star matter with /3-equilibrium 
and charge neutrality as a function of baryon density p obtained from the EOSs listed above 
(Paris+TBF is not shown since it is almost the same as AV18+TBF). For comparison, P for 
the neutron star matter without hyperons obtained from APR EOS ([40]) is also plotted in 
Fig. 2 by the dotted lines. In Fig. 3, the M — pc relationships for corresponding NS models are 
shown. The filled circle on each curve denotes the threshold density of F-mixture. There are 
some features to be remarked from the figure: (i) Different H-EOSs do not show significant 
difference in P up to 2.5 po, and (ii) the y-mixture is delayed from (2-3)po to 4po if there exits 
repulsive three-body force universally for baryons as in the case of TNI2u and TNI3u. Even 
light-mass NSs (M < Mq for TNI2, TNI3 and AV18+TBF and M < 1.2Mq for SCL3A5]) 
have already the F-mixed core. 

In Table 1, we show k, and the threshold densities of hyperon— mixing, pth(A) and pth(S~), 
for each H-EOS. In the same table, we show the maximum-mass M^ax, the radius R and the 
central density pc of the NS obtained from each H-EOS. The values obtained by switching off 
the y-mixing are given in the parentheses for comparison. For the H-EOSs without univer- 
sal three-body repulsion, significant softening due to l^-mixing reduces Mmax, i-e., Mmax 
(without y)=(1.62 - 2.00)-^ M^ax (with y)=(1.08 - 1.26). This clearly contradicts the 
observed mass Mobs = IAAMq for PSR1913+16. On the other hand, H-EOSs with universal 
three-body repulsion (TNI2u, TNI3u), M^ax is recovered nearly to that without Y. 
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Table 1 Properties of various hadronic EOSs with hyperons; TNI2, TNI3, TNI2u, TNI3u [33, 34], Paris+TBF, 
AV18+TBF [36-38] and SCL3AE [39]. k is the nuclear incompressibihty and pth is the threshold density of hyperon- 
mixing with po (=0.17/fm'^) being the normal nuclear density. R and pc denote the radius and central density for the 
maximum mass (Mmax) NS, respectively. The numbers in the parentheses are those without hyperons. *s indicate that 
the numbers are read from the figures in [36]. 



EOS 


TNI2 


TNI3 


TNI2u 


TNI3u 


Paris+TBF 


AV18+TBF 


SCL3AE 


K (MeV) 


250 


300 


250 


300 


281 


192 


211 


Pth(A)/po 


2.95 


2.45 


4.01 


4.01 


2.9* 


2.8* 


2.24 


Pth(S")/po 


2.83 


2.23 


4.06 


4.01 


1.9* 


1.8* 


2.24 




1.08 


1.10 


1.52 


1.83 


1.26 


1.22 


1.36 




(1.62) 


(1.88) 






(2.06) 


(2.00) 


(1.65) 


R(km) 


7.70 


8.28 


8.43 


9.55 


10.46 


10.46 


11.42 




(8.64) 


(9.46) 






(10.50) 


(10.54) 


(10.79) 


Pc/PO 


16.10 


13.90 


11.06 


8.26 


7.35 


7.35 


6.09 




(9.97) 


(8.29) 






(6.47) 


(6.53) 


(6.85) 



The use of several kinds of EOS mentioned above, from different theoretical methods (G- 
matrix, RMF), with various stiffness ranging from k ^ 190 MeV to 300 MeV and with the 
variation of Pth(^) — (2 — 4)po, is expected to cover the present uncertainties of the H-EOSs. 




Fig. 2 Pressure (P) for the y-mixed neutron star matter with /3-equilibrium and charge neutrality as a function 
of the total baryon density p for different types of EOS. Solid red lines: TNI2u (G-matrix approach, universal three- 
body force, hi = 250MeV) and TNI2 (G-matrix approach, three- nucleon force, = 250MeV). Solid blue lines: TNI3u 
(G-matrix approach, universal three-body force, = 300MeV) and TNI3 (G-matrix approach, three-nucleon force, 
K = 300MeV) [33, 34]. Solid green lines: AV18+TBF (G-matrix approach, three-nucleon force, k = 192MeV) [36] and 
SCL3AE (relativistic mean field model with chiral SU(3) symmetry, = 211MeV) [39]. Paris+TEF is not plotted here 
because it is almost the same as AV18+TBF. For comparison, P for the neutron star matter without hyperons obtained 
from APR EOS [40] is also plotted by the dotted lines. 
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Fig. 3 M-pc relationship corresponding to EOSs in Fig. 2 (details of EOS are given in Table Al of Appendix A). 
Colors on each line are the same with those in Fig. 2. The cross symbols denote the points where the NS mass becomes 
maximum, Mmax- Solid black hne denotes M = IAAMq for PSR 1913+16. 



For completeness, numerical values of the pressure P and the energy density 5 as a function 
of the baryon density are tabulated in Table Al in Appendix A. 

3. Quark EOS (Q-EOS) 

The baryon density at the central core of the NSs would be at most lOpo- Although hadrons do 
not keep their identities in such a high density, the chemical potential of the quarks are about 
(400 — 500) MeV which is not high enough for the asymptotic freedom at work. Namely, the 
deconfined quarks inside the NSs, even if they exist, would be strongly interacting. Analogous 
situation at finite temperature has been expected theoretically and is recently confirmed by 
the relativistic heavy-ion collisions at RHIC and LHC; it is now called the strongly interacting 
quark-gluon plasma (sQGP). 

Since lattice QCD to treat the strongly interacting quark matter (sQM) at finite baryon 
density is unfortunately not possible due to the notorious sign problem, we adopt an effective 
theory of QCD, the (2+l)-flavor Nambu— Jona-Lasinio (NJL) model. This model is particu- 
larly useful to take into account the important phenomena such as the partial restoration of 
chiral symmetry at high density [41-44]. 

The model Lagrangian we consider is 

1 ^ 
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where the quark field qi {i = (i, s) has three colors and three flavors with the current quark 
mass rrii. The term proportional to is a U{3)l x U{3)ji symmetric four-fermi interaction 
where are the Gell-Mann matrices with = ^2/3 1. The term proportional to is the 
Kobayashi— Maskawa— 't Hooft (KMT) six-fermi interaction which breaks U{1)a symmetry. 
We consider two types of vector interaction (the second line of Eq.(l)): The term proportional 
to gy{> 0) gives a universal repulsion among different flavors, while the one proportional to 
Gy{> 0) gives flavor-dependent repulsion. 

In the mean-field approximation, the constituent quark masses Mi (i = s) are 
generated dynamically through the NJL interactions {Gs,d), 

Mi = mi- 2G,Gi + 2G^GjGk, (2) 

where di = {qiqi) is the quark condensate in each flavor, and (i, j, k) corresponds to the 
cyclic permutation of u^d and s. The thermodynamic potential is related to the pressure 
as = — TlogZ = — PV, so that we have 

-a,Y, "1 - + 1 "f (3) 

where rii = {qlqi) is the quark number density in each flavor, and Si is the quark propagator, 
which can be written as 

Sr^=^-M,-^yf, t^f^h-';^^""^ (4) 

[fli - G^rii 

where iuj£ = {2£ + l)7rT and fif^ is an effective chemical potential [45]. 

There are six independent parameters in the (2+l)-flavor NJL model; the UV cutoff. A, the 
coupling constants, and ^^(G^), and the quark masses, rriu^d ^s- Five param- 

eters except for g^{G^) have been determined from hadron phenomenology. We consider 
three parameter sets summarized in Table 2; HK (Hatsuda and Kunihiro), RKH (Rehberg, 
Klevansky and Hufner) and LKW (Lutz, Klimt and Weise) [41-44]. 



Table 2 Parameter sets of (2+l)-flavor NJL model [41-44]. 





A(MeV) 


G,A2 


G,A5 


m„,d(MeV) 


m3(MeV) 


HK 


631.4 


3.67 


9.29 


5.5 


135.7 


RKH 


602.3 


3.67 


12.36 


5.5 


140.7 


LKW 


750 


3.64 


8.9 


3.6 


87 



The magnitude of g^{G^) has not been determined well: Recent studies of the PNJL model 
applied to the QCD phase diagram suggest that may be comparable to or larger than G^ 
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Fig. 4 The number fractions {jiu^d.s.el'^tot with 
^T-tot = TT-u + TT-d + ris = 3p) as a function of the baryon 
density p. SoHd red Hne: The fraction of u quark. SoHd 
blue Hne: The fraction of d quark. SoHd green Hne: The 
fraction of s quark. SoHd black line: The fraction of elec- 
tron X 100. Muon does not appear due to the emergence 
of s quarks. 
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Fig. 5 The constituent quark masses {Mi) as a 
function of p. Colors on each line are the same with those 
in Fig. 4. 



[46, 47], so that we change its magnitude in the fohowing range, 

< 1^ < 1.5. (5) 

In §4 and §5, we wih show our results mainly for the HK parameter set with the vector 
interaction of the type. Later, we discuss how the results change in other cases. The 
Q-EOS with strangeness is obtained from the above model under two conditions: (i) the 
charge neutrality among quarks and leptons, i.e. \nu — — \ns — — — and (ii) 
the ^^-equilibrium among quarks and leptons, i.e. fi^ = fJ^s — fJ^u + Me and jie — l^p- 

In Fig. 4, the number fractions {nu^d,s,e/'^tot with ntot = + '^d + = 3p) as a function 
of the baryon density p are plotted. Also, in Fig. 5, the constituent quark masses (Mi) as a 
function of p are plotted. The HK parameter set with the ^^-type interaction are used in both 
figures. The flavor-independent ^^-type interaction leads to a pressure in Eq.(3) depending 
only on /n^^ . Then, the number fractions and the quark masses as a function of p do not 
depend on . 

At low baryon densities below a threshold density pth — 4po, the system is composed of 
only u^d and e with ^ 2nu due to charge neutrality and /3-equilibrium (Fig. 4). In this 
region, the strong interaction among quarks (mainly the G^-term in the NJL model) drives 
the partial restoration of chiral symmetry and hence a rapid decrease of the constituent 
masses M^^d (Fig. 5). Due to the coupling between different flavors through the G^^-term, the 
strange quark mass Mg in the Dirac sea is also affected slightly. 

When the baryon density exceeds pth, the chemical potential of the strange quark ps 
becomes larger than the strange quark mass [ps > Ms)-, so that the system starts to have the 
strangeness degree of freedom. Since the strange quark is negatively charged, the electrons 
start to disappear from the system and the d quark fraction gets decreased at the same time 
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Fig. 6 Pressure (P) as a function of baryon density p in a pure quark matter for the HK parameter set with 
9v l^s ~ ^' ^•^•> "^^^ filled circles denote the onset of the strangeness. 



(Fig.4). We note that the system does not have the muon, because the electron chemical 
potential is smaller than m^=106MeV at all densities. In the high density limit, system 
approaches to the flavor symmetric d, s matter without leptons. Once the 5-quark appears 
in the system, Ms is also suppressed mainly due to the G^-term (Fig. 5). The strangeness 
threshold pth does not depend on as already mentioned, while it depends on the NJL 
parameter sets in Table 2; pth/po — 4.0,3.9 and 3.0 for HK, RKH and LKW, respectively. 

In Fig. 6, we plot the pressure(P) of the strongly interacting quark matter for the HK 
parameter set with different values of the vector coupling {g^ /G^ — 0, 1.0, 1.5 according to 
Eq.(5)). Due to the universal repulsion of the ^^-type vector interaction, the Q-EOS becomes 
stiffer as g^ increases. As mentioned already, the onset density of the strangeness (marked 
by the filled circles) does not depend on g^. We note here that the present Q-EOS has a 
first-order phase transition below 2po for g^ < O.SGg. However, it does not affect the final 
results of the present paper, since such a low density region is dominated by the hadronic 
EOS in our hadron-quark crossover approach to be discussed in §4. 



4. Hadron- Quark crossover 

As discussed in §1, treating the point-like hadron as an independent degree of freedom loses 
its validity as the baryon density approaches to the percolation region. In other words, the 
system cannot be described neither by an extrapolation of the hadronic EOS from the low- 
density side nor by an extrapolation of the quark EOS from the high-density side. Under 
such situation, it does not make much sense to apply the Gibbs criterion of two phases I and 
II, Pi(Tc^fic) — Pu{Tc^fic) since Pi and Pu are not reliable in the transition region. 

Since the first principle QCD calculation at high baryon density is not available and effective 
models at finite baryon density with proper treatment of the confinement phenomena do not 
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exist at present, we will consider a phenomenological "interpolation" between the H-EOS and 
Q-EOS as a first step. Such an interpolation is certainly not unique, but we adopt a simplest 
one^; 

P(p) = PH{p)f-{p) + PQ{p)f+{pl (6) 

/±(p) = l(^l±tanh(^)), (7) 

where Ph and Pq are the pressure in the hadronic matter and that in the quark matter, 
respectively. The window p — V % p + V characterizes the crossover region in which both 
hadrons and quarks are strongly interacting, so that neither pure hadronic EOS nor pure 
quark EOS are reliable. The percolation picture illustrated in Fig.l is best implemented by 
the interpolation in terms of the baryon density p instead of the baryon chemical potential. 

One should not confuse Eq.(7) with the pressure in the mixed phase associated with the 
first-order phase transition in which f± is considered to the volume fraction of each phase. 
In our crossover picture, the system is always uniform and /_ (/+) should be interpreted as 
the degree of reliability of H-EOS (Q-EOS) at given baryon density. 

To calculate the energy density e as a function of p in thermodynamically consistent way, 
we integrate the thermodynamical relation, P = p^d[e/ p)/dp and obtain 

e{p) = eH{p)f-{p)+eQ{p)f+{p) + Ae (8) 
Ae = p [\eHip')-eQip'))^dp' (9) 

J p P 

with g{p) — ^(e^ + e~^)~^ and X — {p — p)/V. Here er (sq) is the energy density obtained 
from H-EOS (Q-EOS). is an extra term which guarantees the thermodynamic consistency. 

In the following, we consider crossover window which satisfies the following physical con- 
ditions: (i) The system is always thermodynamically stable dP/dp > 0, and (ii) the normal 
nuclear matter is well described by the H-EOS so that p — 2T > po is satisfied. Shown in 
Fig. 7 is an example of the interpolation between TNI2u for H-EOS and NJL with = 
for Q-EOS according to Eq.(7). The crossover window is chosen to be (p, F) = (3po,Po) and 
is shown by the shaded area on the horizontal axis. An important lessen one can learn from 
Fig. 7 is that the H-EOS (Q-EOS) is nothing more than the asymptotic form of the "true" 
P{p) around p = (/> = oo). Therefore, naive extrapolation of H-EOS and Q-EOS beyond 
their applicability would miss essential physics. 

In Fig. 8, we plot the interpolated EOS using TNI2u and NJL for different values of in 
a wide range of baryon density. The filled circles denote the onset of strangeness degrees of 
freedom, either hyperons or strange quarks. In the next section, we will solve the structure 
of hybrid stars mainly by using these EOSs. 



^ Similar interpolating function f± has been previously considered at finite temperature in [48] . 
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5. Numerical results and discussions 
5.1. Massive hybrid star with strangeness 

We now solve the following Tolman-Oppenheimer-Volkov (TOV) equation to obtain M-R 
relationship by using the EOSs with and without the hadron-quark crossover: 

= _ (M(r) + 4^Pr^)(s + P)(l-2GM(r)/r)-\ 
dr ^ ^ 

nr 

M{t) = / ^^r''^e{r')dr', (10) 
Jo 

where we have assumed the spherical symmetry with r being the radial distance from the 
center of the star. 

In Fig. 9(a), we show the M-R relationship for various H-EOSs with hyperons whose onset 
is denoted by the filled circles. The crosses denote the points where maximum masses are 
realized: In all cases, M^ax does not reach 2Mq due to the softening of EOS by the hyperon 
mixture. 

In Fig. 9(b), we show the M-R relationship with the EOS interpolated between H-EOS 
and Q-EOS: For the H-EOS, we consider the same EOSs as shown in Fig. 9(a), while for the 
Q-EOS, we adopt the HK-parameter set with = as a typical example. The crossover 
window is fixed to be (p, F) = (Spo^Po)- Cases for different parameters in Q-EOS as well as 
for different window parameters are discussed in the next subsection. 

The red lines in Fig. 9(b) correspond to the cases with TNI2u and TNI2, the blue lines 
correspond to TNI3u and TNI3, and the green lines correspond to SCLSAS and AV18+TBF. 
The onset of strangeness and the maximum mass are denoted by the filled circles and the 
crosses, respectively. Irrespective of the H-EOSs, the interpolated EOS can sustain hybrid 
star with M^ax > SM©: A smooth crossover around p ^ 3po and the stiff Q-EOS due to 
repulsive vector interaction are two fundamental reasons behind this fact. Also, we note 
that the radius of the hybrid star with interpolated EOS is in a range i? = (11 =t l)km for 
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Fig. 9 M — R relationships. (a)M — relationships with various H-EOS including hyperons. Solid red lines: 
TNI2u (universal 3-body force with k = 250MeV) and TNI2 (3-nucleon force with K = 250MeV). Solid blue hues: 
TNI3u (universal 3-body force with k = 300MeV) and TNI3 (3-nucleon force with k = 300MeV) [33, 34]. Solid green 
lines: AV18+TBF (G-matrix approach with hyperons) [36] and SCL3AE (Relativistic mean field model with a chiral 
SU(3) symmetry) [39]. The gray band denotes M = (1.97 ± O.O4)M0 for PSR J1614-2230. The solid black line denotes 
M = 1.44M0 for PSR 1913+16. (b) M-R relationship with the EOS interpolated between H-EOS in (a) and Q-EOS 
with the HK parameter set and = G^, by the window parameters (p, T) = (3po,po)- Colors on each line are the 
same with those in (a). 



0.5 < M/Mq < 2.0, except for the case of SCLSAS. ^ Such a narrow window of R independent 
of the values of M is consistent with the phenomenological constraints on R based on recent 
observations of both transiently accreting and bursting sources [49, 50]. 

In Table 3, we show the maximum mass and the associated central density of the hybrid 
star with the interpolated EOS with g^^G^ and ^^=1.5G^. In all combinations of H-EOS 
and Q-EOS, M^ax exceeds 2M0 with the central density, pc = (5.4 — 6.1)po- 

Let us now turn to the internal structure of the hybrid star, in particular its strangeness con- 
tent. From the location of the filled circles in Fig. 9(b), one finds that the flavor-independent 
universal three-baryon repulsion in TNI2u and TNI3u increases the onset density of the 
strangeness inside the hybrid star. This can be seen more explicitly by plotting the radial 
profile of the hybrid star: The upper panels of Fig. 10 show the p — r relationships for 2Mq 
and 1.44M0 hybrid stars with TNI2 (left) and TNI2u (right). The threshold densities of the 
strangeness given in Table 1 are indicated by the double lines. In our interpolated EOSs, the 



^ The reason why the case with SCLSAS is different from others can be easily seen from Fig. 2: 
The pressure P of SCL3AE is nearly twice as large as that of the other EOSs at p = (1 — 2)po- This 
leads to a larger R of light NSs. 
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Table 3 Mma^/MQ (pc/po) for different choice of H-EOS and difference stiffness of Q-EOS. 



H-EOS 


9v=Gs 




TNI2u 


2.05 (6.1) 


2.17 (5.5) 


TNI2 


2.04 (6.1) 


2.16 (5.9) 


TNI3u 


2.07 (5.9) 


2.18 (5.4) 


TNI3 


2.04 (6.1) 


2.16 (5.5) 


Paris+TBF 


2.06 (6.1) 


2.17 (5.6) 


AV18+TBF 


2.06 (6.1) 


2.17 (5.5) 


SCL3AE 


2.06 (5.9) 


2.17 (5.5) 



TNI2 TNI2U 




Fig. 10 Density-profiles p(r) with r being the distance from the center for the 2.0Mq star (red Hne) and IAAMq 
star (blue line). In the left upper panel TNI2 H-EOS and Q-EOS with = G ^ and the HK parameter set are used 
in the interpolation with the window parameters (p, T) = (Spo^Po), while in the right upper panel TNI2u H-EOS and 
Q-EOS above are used. Double line shows the density above which the strangeness appears. Lower illustrations show 
the internal structure. Only the shaded regions contain strangeness degrees of freedom. 

above stars turn out to have almost the same radius. The lower illustrations of Fig. 10 show 
the cross sections of the corresponding hybrid stars. 

These figures imply that, even if the mass and the radius are the same, the strangeness 
content of the hybrid stars can be quite different. This point is of particular interest for the 
cooling problem of NSs. As is well known, NSs with a y-mixed core undergo an extremely 
rapid cooling due to the efficient //-emission processes called "hyperon direct URCA" (y- 
Durca, e.g., A^p + e~ + z/e, p + e~^A + z/e) and are cooled very rapidly below the 
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detection limit of thermal X-ray. Therefore, for the NSs consisting of pure hadronic com- 
ponents with y, only the very light-mass NSs (M < (1.0 — 1.2) Mq, as in Fig. 3) can escape 
from y-Durca rapid cooling. This means an unlikely situation that all the NSs whose Tg 
are observed should be light-mass stars in spite of the fact that the observed mass distri- 
bution is centered around (1.4 — 1.5)Mq [2]. On the contrary, in the case of the hybrid star 
with = Gg{1.5Gg) under consideration, NSs as heavy as up to 1.9(2.O)M0 can avoid this 
rapid cooling, allowing the T^-observed NSs to be from the light-mass to heavy mass stars 
(M < (1.9 - 2.O)M0, as in Fig. 11). 

It is in order here to comment on the relationship between the maximum mass and the 
nuclear incompressibility hc. From the properties of finite nuclei, the nuclear incompressibility 
hi is estimated to be (240 ± 20)MeV [51]. The interpolated EOSs with TNI2 and TNI2u are 
consistent with this empirical and yet they can reach M^ax > 2Mq. In other words, what 
is important to sustain massive hybrid stars is not the value of the incompressibility, but the 
stiffness of the EOS at and above ^ 3po- 



5.2. Dependence on Q-EOS 




Fig. 11 (a) M — pc relationships with the interpolated EOSs. We adopt the HK-parameter set for the Q-EOS 
with various g^/G^ = 0, 1.0, 1.5. The crossover window are fixed to be (p, F) = (3po,Po)- The cross symbols denote 
the points of Mmax, while the filled circles denote the points beyond which the strangeness appears. The gray band 
denotes M = (1.97 ± O.O4)M0 for PSR J1614-2230. The solid black fine denotes M = I.MMq for PSR 1913+16. (b) 
M-R relationships with the interpolated EOSs. 



To see how the hybrid star structure changes by the stiffness of Q-EOS, we plot M — pc 
relationship for /G^ — 0, 1.0, 1.5 with HK parameter set in Fig. 11(a). We take TNI2u for 
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H-EOS and the same crossover window as in Fig. 9. For comparison, the M — pc relationship 
only with TNI2u is plotted by the dashed line. Fig. 11(b) shows the corresponding M — R 
relations. As anticipated, Mmax increases as increases. In Table 4, we show how Mmax 
and Pc depend on the choice of and the choice of the NJL parameter set. Although 
the parameter dependence is not entirely negligible, the massive hybrid star is possible for 
sufficiently large values of g^ . 

Table 4 The values of Mmax/M© (pc/po) for g^/Gs = 1.0,1.5,2.0 with (p, T) = (3po,po) and TNI2u. The 
parameter sets of the NJL model, HK, RKH and LKW, are given in Table 2. 



Q-EOS 


9v^Gs 




9v=2Gs 


HK 


2.05 (6.1) 


2.17 (5.5) 


2.24 (5.4) 


RKH 


1.99 (6.2) 


2.12 (5.8) 


2.20 (5.4) 


LKW 


1.72 (7.5) 


1.87 (6.7) 


1.97 (6.3) 



Finally, we consider the flavor-dependent vector interaction proportional to given in 
Eq.(l). In the high density limit where c?, s quarks have equal population, (u^u) = (d'^d) = 
(5^5), the g^ interaction and the interaction have the same contribution to the pressure 
in the mean-field approximation if we make the identification, G^ = ^g^. Motivated by this 
relation, we show M^ax and pc for G^/G^ = 1.5, 2.25, 3.0 in Table 5. For the density relevant 
to the core of the hybrid stars, the flavor SU(3) limit is not yet achieved due to the 5-quark 
mass (see Fig. 4). Therefore, the EOS for the flavor-dependent repulsion with G^ — ^9v 
softer than the flavor independent repulsion with g^. This can be seen by comparing the 
corresponding values in Table 5 and those in Table 4. In any case, the massive hybrid star is 
possible for sufficiently large values of G^ . 

Table 5 Mmax/M© {pc/po) for the HK parameter set with the flavor-dependent repulsion Gg. The crossover 
window is (p, F) = (3po,po) and the hadronic EOS is TNI2u. 



Gy — 1.5(j^g 


Gy — 2.25(j^g 


Gy — 3.0G g 


1.87 (6.6) 


1.99 (6.2) 


2.07 (5.8) 



5.3. Dependence on crossover window 

In Table 6, we show Mmax and pc for different choice of the crossover window parameterized 
by p and F. TNI2u and HK parameter set are adopted for H-EOS and Q-EOS, respectively. 
As the crossover window becomes lower and/or wider in baryon density, the interpolated 
EOS becomes stiffer and M^ax becomes larger. To be compatible with the observed massive 
NS with M = (1.97 zb O.O4)M0, the crossover needs to occur in (2 - 4)/>o- 

5.4' Sound velocity of interpolated EOS 

One of the measures to quantify the stiffness of EOS is the sound velocity v^ = ^JdP/de. In 
Fig. 12, we plot Vg for our interpolated EOS with g^/G^ = 0, 1.0, 1.5 as a function of p. The 
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Table 6 Mma.yi/ (pc/po) Under the variation of the parameters, p and F, which characterize the crossover 
window. H-EOS and Q-EOS are obtained from TNI2u and HK parameter set, respectively. Columns without numbers 
are the excluded cases corresponding to p — 2r < po in §4. 



p 


r/po = 1 


r/po = 2 


9v=Gs 




9v=Gs 




3/90 

4po 
5po 
6po 


2.05 (6.1) 
1.89 (7.2) 
1.73 (8.2) 
1.60 (9.6) 


2.17 (5.5) 
1.97 (6.8) 
1.79 (8.0) 
1.64 (9.3) 


1.74 (8.0) 
1.62 (9.2) 


1.80 (7.7) 
1.66 (9.0) 



kinks of i;^ at p ^ 4po are caused by the softening of EOS by the appearance of strangeness. 
The enhancement of of the interpolated EOS relative to the pure hadronic EOS takes 
place just at and above the crossover window. 

5. 5. Stability of hybrid star 

The neutron star is gravitationally stable if the average adiabatic index T satisfies the 
inequality [52]: 



r 4 , GM 
>3+'^- 



(11) 



Here F = din P/dlne is the adiabatic index. Also, XGM/R with A being a numerical constant 
of order unity is a general relativistic correction whose magnitude is much less than 1. Since 



o 
_o 

> 

C 
O 



0.8 - 



0.6 



0.2 



1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 




// Sv=V^^ 






// 3^° ' 

1/ /.-••'' TNI2u 




1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 
P/ P 



5 5.5 6 6.5 7 7.5 







Fig. 12 Sound velocity as a function of baryon density p. Solid lines: Vg — p obtained from the interpolated 
EOS between the H-EOS with TNI2u and the Q-EOS with Qy/Gg = 0, 1.0, 1.5. The crossover window is (2 - 4)po- 
Dotted line : that for pure H-EOS with TNI2u. The filled circles denote the points beyond which strangeness starts to 
appear. 



15/21 



r of our H-EOS is about 2 at all densities and F of our Q-EOS is larger than 4/3 due to the 
constituent quark mass and the repulsive vector interaction, Eq.(ll) is always satisfied and 
our hybrid star is gravitationally stable. 

6. Summary and concluding remarks 

Recent observation of a two-solar mass NS presents a challenging problem how to reconcile 
stiff EOS suggested from the observational side with soft EOS due to hyperon-mixing from 
the theoretical side. In this paper we have studied this problem on the basis of the percolation 
picture from the hadronic matter with hyperons to the quark matter with strange quarks. 
We have constructed an EOS by the interpolation between the H-EOS at lower densities and 
the Q-EOS at higher densities, and found that the hybrid stars could have M^ax ^ 2M0, 
compatible with the observation. This conclusion is in contrast to the conventional EOS for 
hybrid stars derived through the Gibbs construction in which the resultant EOS becomes 
always softer than hadronic EOS and thereby leads to smaller Mmax- 

Our qualitative conclusion is insensitive to the choice of different types of H-EOS and 
different types of vector interaction in Q-EOS, as far as (i) the crossover between the hadronic 
matter and the quark matter proceeds in a relatively low density region, e.g., p = (2 — 4)po, 
and (ii) the quark matter is strongly interacting [g^ /G ^ ^ {1 — 2)). We found that the the 
sound velocity v^^ which increases rapidly in the crossover window for Qy/G^ > 1, can nicely 
characterize the stiffening of the interpolated EOS and associated enhancement of Mmax- By 
increasing further, we can have hybrid stars close to the causality limit, M^ax ~ 2.4M0. 

The idea of rapid stiffening of the EOS starting from 2po opens a possibility that the 
experimental nuclear incompressibility = (240 ± 20)MeV at p ^ po is compatible with the 
existence of massive neutron stars. Also, the idea may well be checked by independent 
laboratory experiments with medium-energy heavy-ion collisions. 

Although the M-R relationship and Mmax are insensitive to the existence of the universal 
three-body repulsion, the onset density of strangeness is rather sensitive to such repulsion. 
If we have three-body repulsion acting universally among baryons, most of the hybrid stars 
with M < (1.9 — 2.0) Mq are free from the extremely efficient hyperon direct-Urca cooling 
process and can avoid contradiction to observations. 

Finally, we remark that the crossover region may contain richer non-perturbative phases 
such as color superconductivity, inhomogeneous structures and so on [1]. How these structures 
as well as the associated cooling processes affect the results of the present paper would be 
an interesting future problem to be examined. 
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AppendixA. EOS tables 

In this Appendix, we show the concrete values of pressure P and energy density s as a func- 
tion of a baryon density x = p/po for H-EOSs, Q-EOSs and interpolated EOSs. 
In Table Al, we show H-EOSs with hyperons; TNI2, TNI2u, TNI3, TNI3u [33, 34], 
AV18+TBF [36] and SCL3AS [39]. 

Shown in Table A2, NJL Q-EOSs with HK parameter set for various vector interactions 
Qy/Gg — 0, 1.0, 1.5 are listed in each row [43]. 

In Table A3, interpolated EOSs between TNI2u H-EOS and NJL Q EOS with HK parameter 
set with Qy/Gg = 0, 1.0, 1.5 for (p, F) = (3po5Po) ^i"^ listed. 
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Table Al Pressure P(MeV/fm3) and energy density £(MeV/fm^) as a function of a baryon density x = p/po 
for various hadronic EOSs with hyperons[33, 34, 36, 39]. 



X 


TNI2 


TNI3 


TNI2u 


TNI3u 


P 


£ 


P 


e 


P 


£ 


P 


£ 


1.0 


2 


162 


3 


162 


2 


162 


3 


162 


1.5 


7 


245 


9 


246 


7 


245 


9 


246 


2.0 


15 


330 


20 


332 


15 


330 


20 


332 


2.5 


26 


418 


32 


422 


26 


418 


38 


422 


3.0 


40 


508 


40 


513 


42 


508 


61 


516 


3.5 


48 


600 


48 


606 


62 


601 


92 


615 


4.0 


57 


693 


58 


700 


87 


697 


131 


718 


4.5 


67 


787 


69 


795 


112 


797 


174 


827 


5.0 


80 


882 


82 


892 


140 


899 


224 


941 


5.5 


93 


979 


96 


990 


173 


1005 


283 


1060 


6.0 


109 


1078 


113 


1090 


211 


1113 


353 


1185 


6.5 


127 


1177 


1 O 1 

131 


iiyi 


255 


1226 


433 


1316 


7.0 


147 


1278 


152 


1293 


304 


1341 


525 


1454 


7.5 


169 


1381 


175 


1397 


360 


1461 


629 


1599 


8.0 


193 


1485 


200 


1503 


422 


1584 


746 


1751 


8.5 


219 


1590 


227 


1610 


492 


1712 


876 


1911 


9.0 


248 


1698 


256 


1719 


568 


1843 


1020 


2079 


9.5 


278 


1807 


287 


1829 


651 


1979 


1178 


2256 


10.0 


311 


1917 


321 


1942 


743 


2120 


1352 


2441 



X 


Paris+TBF 


AV18+TBF 


SCL3AI] 


P 


£ 


P 


£ 


X 


P 


£ 


0.471 


0.370 


71.2 


0.432 


71.2 


1.02 


5.85 


167 


0.941 


2.41 


148 


2.59 


148 


1.51 


14.6 


251 


1.18 


4.63 


187 


4.94 


187 


2.04 


28.9 


346 


1.76 


15.5 


287 


15.5 


288 


2.51 


41.5 


434 


2.35 


31.2 


404 


29.0 


399 


3.02 


53.5 


531 


2.94 


45.7 


503 


44.2 


505 


3.54 


67.8 


635 


3.53 


62.3 


617 


59.9 


623 


4.07 


84.0 


740 


4.12 


79.0 


735 


75.3 


729 


4.57 


101 


841 


4.71 


99.9 


853 


94.4 


858 


5.01 


118 


933 


5.29 


117 


965 


112 


970 


5.62 


144 


1063 


5.88 


145 


1128 


139 


1122 


6.02 


162 


1150 


6.47 


168 


1240 


159 


1223 


6.60 


191 


1278 


7.06 


188 


1369 


181 


1341 


7.07 


216 


1384 


7.65 


213 


1470 


205 


1459 


7.58 


244 


1499 


8.24 


242 


1599 


239 


1616 


8.12 


276 


1625 


8.82 


279 


1745 


270 


1745 


8.50 


300 


1715 


9.41 


307 


1874 


302 


1879 


9.11 


339 


1860 


10.0 


347 


2031 


328 


1986 


9.54 


368 


1965 
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Table A2 Pressure P(MeV/fm^) and energy density £(MeV/fm^) as a function of a baryon density x = p/po 
for NJL Q-EOSs with HK parameter set for = Gg [43]. 



X 


9v/Gs 


— 


9v/Gs = 1 


9v/Gs 


= 1.5 


r 


£ 


P 


£ 


r 


£ 


1.0 


-0.7633 


179.2 


8.390 


188.3 


12.97 


192.9 


1.5 


-1.397 


268.1 


19.20 


288.7 


29.49 


299.0 


2.0 


6.721 


357.8 


43.33 


394.5 


61.64 


412.8 


2.5 


28.71 


451.3 


85.91 


508.5 


114.5 


537.1 


3.0 


58.56 


550.0 


140.9 


632.4 


182.1 


673.6 


3.5 


91.96 


654.0 


204.1 


766.2 


260.1 


822.2 


4.0 


127.5 


762.8 


274.0 


909.3 


347.2 


982.5 


4.5 


158.2 


876.1 


343.6 


1061 


436.2 


1154 


5.0 


182.8 


992.3 


411.6 


1221 


526.0 


1335 


5.5 


203.0 


1111 


479.9 


1388 


618.3 


1526 


6.0 


221.0 


1231 


550.5 


1561 


715.2 


1725 


6.5 


238.9 


1353 


625.6 


1740 


819.0 


1933 


7.0 


258.3 


1476 


706.8 


1925 


931.0 


2149 


7.5 


279.8 


1601 


794.7 


2115 


1052 


2373 


8.0 


303.7 


1727 


889.5 


2313 


1182 


2605 


8.5 


330.0 


1854 


991.3 


2516 


1322 


2846 


9.0 


358.4 


1984 


1100 


2725 


1471 


3096 


9.5 


389.9 


2114 


1215 


2941 


1628 


3354 


10 


421.3 


2247 


1337 


3163 


1794 


3620 
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Table A3 Pressure P(MeV/fm3) and energy density £(MeV/fm^) as a function of a baryon density x = p/po 
for interpolated EOSs between TNI2u H-EOS and NJL Q EOS with HK parameter set with g^/G^ = 0, 1.0, 1.5 for 
(p, r) = {3po,po) case. For the practical use, in a hadronic phase, we fit pressure as a function ap^ by using the 
least-squares method. 



X 


TNI2 


TNI2u 


gJGs = 


9v/Gs = 1 


9v/Gs 


= 1.5 


= 


9vlGs = 1 


gJGs 


= 1.5 


p 


e 


P 


e 


P 


e 


p 


e 


P 


e 


p 


e 


1.0 


2.480 


168.9 


2.645 


179.5 


2.727 


184.8 


2.480 


168.8 


2.645 


179.4 


2.727 


184.8 


1.5 


6.737 


255.4 


7.714 


271.5 


8.202 


279.6 


6.737 


255.3 


7.714 


271.4 


8.202 


279.5 


2.0 


13.90 


343.7 


18.26 


366.0 


20.45 


377.1 


13.90 


343.6 


18.26 


365.9 


20.45 


377.0 


2.5 


26.93 


434.4 


42.31 


464.4 


50.00 


479.4 


26.93 


434.3 


42.31 


464.3 


50.00 


479.3 


3.0 


49.03 


528.1 


90.22 


569.3 


110.8 


589.8 


50.19 


528.5 


91.38 


569.7 


112.0 


590.3 


3.5 


80.26 


629.2 


162.2 


687.2 


203.2 


716.2 


83.90 


627.5 


165.8 


685.5 


206.8 


714.5 


4.0 


119.2 


733.0 


248.2 


814.2 


312.7 


854.8 


122.7 


731.6 


251.7 


812.8 


316.2 


853.4 


4.5 


153.9 


841.7 


330.4 


952.0 


418.7 


1007 


156.0 


841.7 


332.6 


952.0 


420.9 


1007 


5.0 


181.0 


953.8 


405.7 


1099 


518.0 


1171 


182.1 


954.0 


406.8 


1099 


519.1 


1171 


5.5 


202.3 


1068 


477.3 


1252 


614.8 


1344 


202.8 


1069 


477.8 


1253 


615.4 


1345 


o.O 


220.7 


1185 


549.4 


1413 


713.7 


1527 


220.9 


1185 


549.6 


1413 


714.0 


1527 


6.5 


238.8 


1303 


625.1 


1579 


818.3 


1718 


238.9 


1303 


625.3 


1580 


818.4 


1718 


7.0 


258.2 


1422 


706.6 


1752 


930.7 


1917 


258.3 


1422 


706.6 


1752 


930.8 


1918 


7.5 


279.8 


1543 


794.6 


1931 


1052 


2125 


279.8 


1543 


794.6 


1931 


1052 


2125 


8.0 


303.7 


1665 


889.5 


2115 


1182 


2341 


303.7 


1665 


889.5 


2116 


1182 


2341 


8.5 


330.0 


1789 


991.3 


2306 


1322 


2565 


330.0 


1789 


991.3 


2307 


1322 


2566 


9.0 


358.4 


1914 


1100 


2503 


1471 


2798 


358.4 


1915 


1100 


2504 


1471 


2799 


9.5 


389.0 


2041 


1215 


2707 


1628 


3039 


388.9 


2042 


1215 


2707 


1628 


3040 


10.0 


421.3 


2170 


1337 


2916 


1794 


3289 


421.3 


2170 


1337 


2917 


1794 


3290 
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